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Competition between electron and phonon excitations in the scattering of nitrogen 
atoms and molecules off tungsten and silver surfaces 
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We investigate the role played by electron-hole pair and phonon excitations in the interaction 
of reactive gas molecules and atoms with metal surfaces. We present a theoretical framework 
that allows us to evaluate within a full-dimensional dynamics the combined contribution of both 
excitation mechanisms while the gas particle-surface interaction is described by an ab-initio potential 
energy surface. The model is applied to study energy dissipation in the scattering of N2 on W(110) 
and N on Ag(lll). Our results show that phonon excitation is the dominant energy loss channel 
whereas electron- hole pair excitations represent a minor contribution. We substantiate that, even 
when the energy dissipated is quantitatively significant, important aspects of the scattering dynamics 
are well captured by the adiabatic approximation. 
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In the last years, with the development of ab-initio 
calculations based on density functional theory (DFT), 
unprecedented accuracy has been achieved in describing 
the interaction of reactive thermal and hyperthermal gas 
molecules and atoms with metal surfaces. In most ad- 
vanced simulations, molecular dynamics are performed 
in ab-initio (DFT) ground state multidimensional poten- 
tial energy surfaces (PES) [l|. This refined scheme relies 
on the validity of the adiabatic approximation. By adia- 
batic, we refer to a process which neglects both electronic 
excitations and lattice phonons excitations. Neverthe- 
less, the role played by these nonadiabatic effects is now 
under close scrutiny. When the experimental results and 
the adiabatic results are at variance, quite often contro- 
versy arises about up to what extent the differences are 
due to the neglect of energy loss channels [2| , to the limi- 
tations of a reduced dimensionality approximation |3f , or 
to the inherent limitations of DFT 

The increasing number of gas-surface experiments re- 
porting evidence of electronic excitations [7|, [8| , has mo- 
tivated the development of different models to include 
this dissipation channel [9l-ll2J. Among them, the local 
density friction coefficient approximation (LDFA) [10( of- 
fers a very good compromise between the accuracy of the 
results and the simplicity of its implementation [13j . Dif- 
ferent approaches have been also developed to treat the 



energy exchange with the lattice [144171 ]. Scmiclassical 
approximations of the phonon excitations, in which the 
gas-surface interaction is limited to simplified model po- 
tentials 17( , have been very successful in understanding 



the scattering of nonreactive rare gas atoms with sur- 
faces [10]. However, they are less accurate when applied 
to reactive species [l9(. In the latter case, the multi- 
dimensional PES is typically so corrugated and intri- 
cate that a more realistic treatment of the interaction 



is required. The generalized Langevin oscillator model 
(GLO ) |1 4l| shows to be a sound alternative in this re- 
spect [201 ] . 

Still, the challenge in gas-surface dynamics is to pro- 
vide a theoretical framework that, keeping the accuracy 
of a multidimensional ab-initio PES for the gas-metal 
interaction, incorporates into the dynamics energy ex- 
change with both lattice vibrations and electronic excita- 
tions [2l|, |22| . In this Letter we accomplish this objective 
by combining the GLO for phonon excitations and the 
LDFA for electronic excitations. The inclusion of both 
effects will allow us to address such fundamental ques- 
tions as (i) what is the relative importance of phonon 
and electron- hole (e-h) pair excitations as energy dissi- 
pation channels, in particular, (ii) is there any coupling 
between them or is their contribution just additive, and 
(iii) to what extent does the adiabatic calculation capture 
the basic physics of the dynamics and provide accurate 
results. 

To answer these questions, we have selected two sys- 
tems for which accurate energy loss measurements exist: 
the rotationally inelastic scattering of N2 on W(110) 12311 
and the scattering of hyperthermal N on Ag(lll) [241 ] . 
Molecular and atomic nitrogen are relatively heavy pro- 
jectiles for which energy exchange to the lattice can be 
important. The use of N also implies that electronic ex- 
citations are analyzed under the most favored conditions 
of open-shell reactive species. 

The way in which we incorporate electronic and 
phonon excitations into the multidimensional classical 
trajectory simulations is as follows 25(. As in the GLO, 
the motion of surface atoms is represented by a three- 
dimensional (3D) harmonic oscillator. The latter is cou- 
pled to a 3D ghost oscillator that is subjected to frictional 
and random forces accounting for energy dissipation and 
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FIG. 1. (Color online) Total final energy normalized to the 
incident energy (left ordinate) and energy loss (right ordi- 
nate) of N2 scattered off W(110) vs the exit rotational en- 
ergy. Results for normal incidence and detection angles and 
two incidence energies Ei. In triangles, the experimental data 
of [23|. Our simulations are represented by red filled circles 
(LDFA), blue open squares (GLO), and green filled squares 
(full nonadiabatic calculation). T s — 1200 K. 



thermal fluctuations [20[ • The effect of electronic excita- 
tions is added by introducing for each of the impinging 
gas atoms a separate friction force proportional to its ve- 
locity. The electronic friction coefficient is calculated at 
each point of the trajectory as that of the atom moving 
in a homogeneous free electron gas with electronic den- 
sity equal to that of the surface at this point [10( . Very 
recently, it has been shown that this model constitutes 
an efficient and sufficiently accurate tool to incorporate 
electronic excitations within multidimensional molecular 
dynamics [26j |. 

When a rotationally cold N2 beam is scattered off 
W(110), the molecules lose around 30% of their inci- 
dence energy for low exit rotational states, but the en- 
ergy loss is lower for those scattered at high exit ro- 
tational states [23[. To understand these findings, we 
have performed classical molecular dynamics simulations 
with different approximations: (i) using the adiabatic ap- 
proximation, i.e. neglecting e-h pair and phonon excita- 
tions, (ii) including only e-h pair excitations (LDFA), 
(iii) including only energy exchange with the lattice 
(GLO) and (iv) including both energy dissipation chan- 
nels (full nonadiabatic calculation with respect to both 
electrons and phonons). In all cases we use the 6D ab- 
initio N 2 /W(110) PES of IH [13, where the PW91 ex- 
change correlation functional was used [29f . A minimum 



of 30000 trajectories is calculated using a conventional 



Monte Carlo sampling of all possible initial conditions. 
In (iii) and (iv), the parallel and perpendicular surface 
oscillator frequencies are 19 and 16 meV respectively [30J, 
and the friction coefficients of the ghost oscillators are 
obtained from the Debye frequency as proposed in [14|. 

In Fig. Q] we compare the results of our simulations 
with the experiments. The figure shows the fraction of 
energy retained and the energy loss as a function of the 
exit rotational energy for the scattered molecules, at nor- 
mal incidence and detection angles. When only e-h pair 
excitations are included, contrarily to the experimental 
observations, the energy loss, which is marginal, does not 
depend on the exit rotational energy. It only depends on 
the total energy no matter how it is distributed among 
the different degrees of freedom. The inclusion of phonon 
excitations changes the picture completely. We recover 
the experimental observation that more energy is lost at 
low exit rotational states. Note also that when phonon 
excitations are included the results with and without e- 
h pair excitations are indistinguishable. The differences 
are, within the statistical errors, of the order of the con- 
tribution of e-h pair excitations alone. This reflects the 
predominant role of phonon excitations in this kind of 
experiment. 

The importance of including energy exchange with the 
lattice can be rationalized by noticing that at normal inci- 
dence and detection, though corrugation and anisotropy 
of the PES may complicate the picture, backscattering 
conditions must prevail. As a consequence, large mo- 
mentum transfer takes place from the projectile to the 
lattice in the direction normal to the surface. This im- 
plies a comparatively larger probability for translational 
energy transfer. As a result, molecules that are rotation- 
ally excited at the expense of their translational energy 
in the scattering with the surface, are more inefficient 
transferring energy to the lattice. In Ref. [15j, a sim- 
plified kinematic model has been proposed to illustrate 
this effect. Note that at the highest rotational states, 
within the experimental error bars, it cannot be decided 
whether the molecules overcome minor energy losses or 
energy gains. In fact, for Ei=0.5 eV our calculations 
give a small energy gain consistent with the experimen- 
tal error bars. In this case, the efficient conversion of 
translational energy into rotational energy implies a big 
reduction of the former, and, for the high temperatures 
under consideration (T s = 1200 K), energy transfer from 
the lattice to the projectile is favored. The e-h pair exci- 
tation mechanism cannot give rise to this kind of behavior 
due to the large mismatch between the projectile and the 
electron masses. 

Interestingly, the measured rotational state population 
distributions of the scattered N2 are already well re- 
produced within the adiabatic approximation, as shown 
in |6j. In Fig. [2] we show that including energy exchange 
with lattice and e-h pair excitations the width and shape 
of the distributions are not altered significantly. In other 
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FIG. 2. (Color online) Final rotational state population distri- 
butions for N2 scattered from W(110) under the same exper- 
imental conditions of Fig. [l] The experimental data of [23J 
(filled triangles) are compared with the adiabatic [0] (open 
circles) and the full nonadiabatic (green filled squares) calcu- 
lations. 



words, the conversion from translational to rotational 
energy in this kind of experiment is basically an adia- 
batic process, and the result of the scattering on a highly 
anisotropic and corrugated six-dimensional PES. Hence, 
even when the molecules lose a significant amount of en- 
ergy, the adiabatic approximation is still valid to describe 
different aspects of the molecule-surface interaction, pro- 
vided an accurate full-dimensional PES is used. The pic- 
ture that emerges from our results is the following: in the 
scattering of the molecules with the anisotropic and cor- 
rugated PES, energy is transferred adiabatically from the 
translational to the rotational degrees of motion. Since 
translational energy is more efficiently transferred to the 
lattice than rotational energy, molecules that remain in 
the lowest rotational states (with higher translational en- 
ergy) are the ones that lose more energy. 

Next, we use our model to analyze the energy loss ob- 
served in the scattering of N atoms with energies of some 
cVs on the Ag(lll) surface [24}. These measurements 
represent an excellent benchmark to explore the accuracy 
of our model for atoms and higher projectile energies, for 
which e-h pair excitations are expected to be more rele- 
vant. Figure [3]reproduces the experimental results for an 
effusive beam of N scattered from Ag(lll) at T s = 500 K 
and an incidence angle Oi = 60°. The beam has an av- 
erage energy of 4.3 eV and a FWHM of - 5.0 eV. The 
figure depicts the ratio between the average final < Ef > 
and initial < Ei > energies as a function of the in-plane 
scattering angle O t (see inset). The experiments show 
a decrease of the final average energy as the scattering 
angle increases. A remarkable feature is the dramatic in- 
crease of the energy ratio at t < 60°, i.e., for grazing 
outgoing angles. In fact, the final average energy is larger 
than the initial one in this angular range. 

We have performed classical dynamics calculations, us- 
ing the ab-initio 3D N/Ag(lll) PES of [H|. The results 
of our different simulations for a monoenergetic beam 
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FIG. 3. (Color online) Ratio of final to initial average en- 
ergy vs the total scattering angle [Q t = 180° — (0i + Of)] for 
N/Ag(lll) at T s = 500 K and 6* = 60°. The experimental 
data of [2J] (open and filled triangles correspond to different 
runs of the same experiment) are compared with our simu- 
lations for a monoenergetic (left panel) and for an effusive 
beam (right panel): adiabatic results (open circles), LDFA 
(red filled circles), GLO (blue open squares), and full nonadi- 
abatic (green filled squares). 



with Ei = 4.3 eV and 0j = 60° arc represented in the 
left panel of Fig. [3J For each kind of simulation, we cal- 
culate 300000 trajectories to assure good statistics. The 
parallel and perpendicular surface oscillator frequencies 
used in the GLO and full nonadiabatic calculations are 14 
and 9 meV respectively [32|. The adiabatic results that 
correspond to < Ef > / < Ei >= 1 are not shown in 
the figure. Despite the large incident energy, we observe 
that e-h pair excitations (LDFA calculations) produce 
marginal energy losses. In addition, the electronic energy 
loss is roughly independent of the scattering angle. When 
energy exchange with the lattice is included (GLO and 
full nonadiabatic calculations), the experimental results 
for large scattering angles are well reproduced. Again, 
differences between the results of these two models are 
of the order of the contribution of c-h pair excitations 
alone. However, none of the energy dissipation channels 
is able to explain the energy ratios larger than 1 that are 
measured for small scattering angles. We find that this 
behavior is indeed due to the effusive beam itself. 

The results of our calculations mimicking the exper- 
imental effusive beam are shown in the right panel of 
Fig. El In this case, we use 600000 trajectories. For large 
scattering angles (0 t > 80°), the results do not signifi- 
cantly differ from those obtained with the monoenergetic 
beam. However, at small scattering angles, all the effu- 
sive beam simulations give final to initial energy ratios 
greater than one. Quantitative agreement with the ex- 
periments is obtained only when energy exchange with 
the lattice is allowed, no matter whether e-h pair exci- 
tations are included or not. The small discrepancy ob- 
served at t ~ 60° might be due to the composition of 
the N experimental beam containing not only the ground 
state N( 4 5) but also electronically excited states N( 2 £>), 



N( 2 P), assumed to be mainly reflected around the spec- 
ular position (Q t = 60°) 0. 

A capture mechanism has been proposed to explain the 
small-0 t behavior |24j |. This capture mechanism assumes 
that the lower energy atoms with grazing exit trajecto- 
ries are trapped close to the Ag surface, resulting in the 
increase of the average final energy. An analysis of our 
simulations shows that the small-Gt behavior is due to 
the fact that the higher energy atoms (Ei > 7 eV) of the 
effusive beam are preferentially reflected at low scatter- 
ing angles. As a result, the average energy of the atoms 
scattered at small (large) 9 t is higher (lower) than the av- 
erage energy of the incident beam. We stress that this is 
a pure trajectory effect that is not related to the presence 
or not of inelastic channels: our adiabatic simulations for 
the effusive beam already reproduce the correct depen- 
dence of < Ef > j < Ei > on 9 t . This clearly indicates 
that a realistic PES is essential in order to reproduce in 
quantitative terms this nontrivial dynamic effect. 

In summary, we have developed a theoretical frame- 
work that allows the inclusion of both phonon and e-h 
pair excitations in a full-dimensional dynamics, keeping 
the accuracy of an ab-initio PES to treat the particle- 
surface interaction. Though the model is not applica- 
ble to systems in which nonadiabatic effects are due to 
crossing of potential energy surfaces with possible charge 
transfer, it can provide a proper description for a great 
variety of systems. Concerning the three questions posed 
above, our analysis shows that for two representative sys- 
tems: (i) phonon excitation is the main energy dissipa- 
tion channel even for the hyperthermal N beams, (ii) the 
contribution of phonon and e-h pair excitations to the 
total energy loss is additive, and, finally, (iii) even when 
energy loss processes are significant, important aspects 
of the scattering dynamics are already captured by the 
adiabatic calculation: the conversion of translational to 
rotational energy upon scattering with the surface in the 
Na/W(110) case and the decrease in the number of low 
energy atoms that are scattered off the surface with small 
grazing exit angles in the N/Ag(lll) case. The correct 
description of these features requires a full-dimensional 
dynamics calculation on top of an accurate PES. In this 
respect, the theoretical framework proposed here is an 
excellent choice to efficiently include all the basic ingre- 
dients in gas-surface simulations. 
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